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Abstract 

A criterion, proposed by the present authors, is used to derive numerical 
results for the stability of a large bipolaron embedded in a polaron gas. The 
main conclusion is that an isolated metastable bipolaron can be stabilized by 
the polaron gas surroundings because of the Fermi statistics of polarons. On 
the other hand, it is found that the exchange interaction tends to destabilize 
the bipolaron. The study is performed both for bulk (3D) materials and for 
thin (2D) films within the Hartree-Fock approximation. The bipolaron is 
described by an extension of the Feynman polaron model. 
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Two electrons (or holes) in a polar crystal interact with the phonon field which may 
lead to their binding into a composite quasi-particle, a bipolaron. The binding into a bipo- 
laron results if the attraction between the electrons, due to the virtual phonon exchange, is 
sufficiently strong to overcome the electron-electron repulsion. 

In this letter we deal with the so called large (bi)polarons (e.g. Refs. for more 

references see the review articles H|7|)- The singlet large bipolaron is characterized by two 
dimensionless coupling constants. The first one, the Frohlich coupling constant 

a = J-S(2--L)J™^ (i) 

nuJ LO V2 V^oo EqJ V h 

is a measure of the strength of the electron-phonon interaction while the second one 

D= 1 *_MW£ 
hLJ LO Eoo V h 

is the Coulomb potential coupling constant which governs the strength of the direct electron- 
electron repulsion. In ([[J) and (0) e and are the static and the high-frequency dielectric 
constants, is the electron band mass and u LO is the frequency of the longitudinal optical 
(LO) phonons. Introducing the ratio of the dielectric constants r] = e^/so, one obtains 
the following relation between the Coulomb and the electron-phonon coupling constants 
U = y2a/(l — As r] > 0, only values U > \/2a have a physical meaning. In what 
follows use is also made of the parameter 

u = _ = . 3 

The physical relevant region is then defined by the inequality u > 1. 

To find the stability region of the bipolaron in the (a,w)-plane one usually exploits the 
following condition 

Ebip < 2E po i, (4) 

where E po i and Eu P denote the ground state energies of the polaron and the bipolaron at 
rest, respectively. This inequality implies that the decay of the bipolaron into two polarons 



is not energetically advantageous. With the bipolaron turns out to be stable if a is larger 
than some critical coupling constant a c . In three dimensions (3D) the following estimates 
have been reported for the critical value of the electron-phonon coupling constant: a c = 7.3 
in Ref. a c — 6.8 in Ref. [[| and a c = 6 in Ref. 0. Even if a exceeds the critical value a c , 
the possibility of large bipolaron formation depends on the strength of the direct Coulomb 
repulsion which tends to prevent electrons to congregate into a cluster. Namely, Eu v obeys 
(|) at r] < r] c (a) (or, equivalently, u < u c (a)). For example, for large coupling, a = 9, it has 
been found that rj c « 0.056 {u c « 1.059) in Ref. § and r] c « 0.037 (u c « 1.038) in Ref. §. 
(The quoted results for a c and r] c , u c leave little hope that the 3D singlet large bipolaron 
can survive even in special kinds of materials (such as strongly ionic crystals)). 

In reality, the inequality (^) is the stability criterion for an isolated bipolaron; it is not 
applicable for a bipolaron interacting with a system of charge carriers. In Ref. || the present 
authors have proposed a stability criterion for a 3D bipolaron embedded in a polaron gas. 
The basic idea is the following. Polarons are fermions and obey the Pauli principle. Therefore 
the two polarons in a final state after the bipolaron decay must have their momenta outside 
the Fermi-surface. Consequently, the total kinetic energy of two such final state polarons 
cannot be less than 2{p F /2m po {) where m po i is the polaron effective mass. We must add this 
term to the r.h.s of the criterion @. This evidently makes the bipolaron decay less probable. 
That is, the bipolaron is possibly stabilized because of the Fermi statistics of the final state 
polarons. Furthermore, the polaron-polaron and the bipolaron-polaron interactions should 
be taken into account. 

In the present letter we restrict ourselves to the Hartree-Fock approximation and perform 
the calculations for an arbitrary number D of dimensions to obtain the following formulae 
for the 2D bipolaron. The .D- dimensional Fermi-momentum is given by the expression 

l/D 
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where n = N/V is the polaron concentration in a D-dimensional box. In the Hartree-Fock 
approximation the mean kinetic energy per particle (without polaron effects) reads as follows 



W(n) = (-?—) 
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A term to be added to the r.h.s of the criterion (J^j is the difference between the kinetic 
energies of N + 2 polarons (after the bipolaron decay) and that of iV polarons (the bipolaron 
is assumed to be at rest): 
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Thus, ^kin is twice the polaron kinetic energy on a Fermi-surface. 

In a similar way we treat the potential energy of the polaron gas. At large distances 
t ^ r poi (where r po i is the polaron radius) the polaron-polaron and the bipolaron-polaron 
interactions can be approximated by Coulomb potentials screened by the static dielectric 
constant: $ pp ~ e 2 je^r and $b p ~ 2e 2 /eor. In the Hartree-Fock approximation the exchange 
energy per particle is then approximated as follows 
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The integral in (||) can be evaluated 

r d D kd D q 

\k\,\q\<PF 

and we obtain from Eq. (|8|) 
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The exchange energy to be added to the r.h.s of the criterion (ffl) takes therefore the form 
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Finally, in the Hartree-Fock approximation, we arrive at the following in-medium crite- 
rion for the bipolaron stability 

Pp c 2 4 pp 

E bip - 2E po i < — — — ■ — . (12) 

m po i e n(D - 1) h 

The r.h.s. of ( |T2"D is equal to zero for p F = and for p F = p d , where p d is given by the 
relation 

Tripoli 4 /,„•, 

Pd = -f^^D^Yy (13) 

In the interval < pp < pd, that is at small polaron concentrations, the r.h.s. of the criterion 
(|i~2|) is negative which allows the bipolaron to decay even if it is stable when isolated. In 
this case the exchange energy dominates and a "stable" bipolaron (in the sense of the 
criterion (^)) is destabilized by the polaron gas. In the interval Pf > Pd the effect of the 
Fermi statistics dominates and the r.h.s. of ([12]) is positive thus additionally stabilizing the 
bipolaron in comparison with what follows from (|4]). 
At D = 3 we arrive at the stability criterion || 

2E pol < — (3vr 2 n) 2/3 - 2 * I*-) 1 ". (14) 
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In dimensionless units the criterion (0) takes the form 

E Uv - 2E vo i m b ( o n \ 2 / 3 n , ( 3 n \ 1/3 , s 

v — < — - 3tt 2 — -2\/2 a . 15 

huj^o nipoi \ n J \tc n J 

Here Uq denotes the natural unit for the polaron concentration 

no=(—) • (16) 

Its typical value can be estimated with the data of Ref. 0. For instance, for RbCl irtb = 
0.432m e (m e is the electron mass) and uj h0 = 3.4 ■ 10 13 s -1 , so n ~ 4.5 ■ 10 19 cm -3 . 
For quasi-particles confined to a two-dimensional (D = 2) layer (O) takes the form 



Ebip - 2E pol <—2na-4- J—, (17) 
mpoi e V 7T 



where we used the notation a for the surface concentration of polarons. Eq. ( ]T7|) can be 
written as follows 
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We used the criteria ( fHf) , ([15|) in Ref. to study the possibility of stabilizing the 
bipolaron, restricting ourselves to the limiting case u = 1 (that is, omitting the exchange 
energy which diminishes the bipolaron stability region). The goal of the present letter is 
to study both stabilization and destabilization effects in more detail using the criteria flHID, 
(|T%|). To calculate the bipolaron and the polaron ground state energies use is made of the 
bipolaron model of Ref. || which is a generalization of the Feynman model for the single 
polaron. As for any Gaussian approximation, the results for the polaron ground state energy 
and the effective mass in different dimensions are linked to each other by a scaling relation 



introduced in Ref. 10 
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A similar scaling relation is valid for the bipolaron energy || (the Coulomb coupling constant 
is not scaled). At the "destabilizing" concentration corresponding to the value pd of Eq. 
(|T3|), the r.h.s. of the stability criterion equals zero. The bipolaron energy then equals twice 
the polaron energy. This occurs at the boundary of the bipolaron stability region found with 
the old criterion (^). In 3D we obtain rid from ([15]) 



and in 2D from (IT 
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To find a simple estimate for nd we may use the interpolation formula derived in Ref. [T2 
for the bipolaron model used here in 3D: 
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u c {a) = 1.08525 V-y^ ( 22 ) 
or — 7.969 

(note that w of the present letter differs by a factor y/2 from w of the paper |~2]]). Inserting 
Eq. (|22|) into Eq. (^) and using the value of the polaron mass calculated in the scope of 
the Feynman model, we may estimate n d . For instance, at a = 7 we obtain m po i/mb = 14.4 
and n d = 0.01n . The strong dependence of n d on a is mostly due to the strong dependence 
of rrip i on it. That is, the smaller is a the less is the bipolaron destabilization effect. 

The interpolation formula for u c (a) in 2D can be found with the same scaling law (|19|): 

uf D \a) = u(? D \3a/4n). It follows then for the 2D case 

a 2 - 1.968 . , 

«.(«) = 1.08525 ^- Tis . (23) 

Because of the scaling relations the value a = 4x9/37r = 3.82 in 2D produces the same 
polaronic effects as a = 9 in 3D. For this a the polaron mass takes the value m po i/mb = 62.7 
and Eq. ( |21~D gives the surface concentration o d ~ 50o"o. The numerical solution leads at 
the value 0^ ~ 43<t . For a = 2.97 in 2D, which corresponds to a = 7 in 3D, we obtain 
a d w 0.04a . 

This means that the bipolaron which is stable when isolated can decay into two polarons 
for concentrations n < n d [a < a d ). In the opposite case n > n d (a > a d ) our in- medium 
stability criterion allows the bipolaron to exist even if its energy exceeds twice the polaron 
energy. Such a bipolaron being isolated would decay, therefore we call it a metastable 
bipolaron. For a given value of a a metastable bipolaron is seen in the range of the Coulomb 
coupling constant u c (a) < u < u max (a). Only a state of two free polarons was seen at 
u > u max (a). 

In Fig. p] the bipolaron ground state energy is shown as a function of the Coulomb 
coupling constant u for a = 9 in 3D. Subsequently, the same curve is obtained in 2D for 
a = 4 x 9/(37r) = 3.82 (cf. the right y-axis). In this case the isolated bipolaron is stable 
for 1 < u < u c = 1.038. In the region u c < u < u max = 1.146 the bipolaron energy exceeds 
twice the polaron energy but this metastable state can be stabilized by the polaron gas. At 
u > u max our numerical program has found only two separate polarons minimum. 



In Fig. |2| regions are shown where polarons, bipolarons and metastable bipolaron states 
exist. The top x-axis corresponds to the 2D case. The metastable bipolaron is seen to occur 
for a > a m i n ~ 6.5 (in 3D) which is slightly smaller than the critical value a c ~ 6.8. Note 
also that a metastable bipolaron was also reported in Ref. [ |TTf . 

It follows from our stability criterion that the metastable bipolaron is stabilized by the 
polaron environment for concentrations n > n^. The stabilizing effect is evident for a 
close to the critical value (we choose a = 7 for the 3D case in Fig. §): the limiting value 
of the Coulomb coupling constant increases with the polaron concentration n. Bipolarons 
can exist for the parameter values below the plotted curve. The stability region of the 
isolated bipolaron is determined by u < u c . At n = nu m pa 2.4n the constant u reaches 
its maximal value u max . Note that typical concentrations for stability of bipolarons ~ 
10 18 -7- 10 19 cm -3 are close to charge-carrier concentrations of heavily doped degenerate polar 
semiconductors pT^— p~5|| . 

The similar plot for a = 3.82 in 2D is presented in Fig. |j. Both the stabilization (at 
a > o d pa 43cr ) and the destabilization (at a < a a) effects are seen clearly. 

To conclude, we have formulated an in-medium stability criterion for bipolarons which 
takes into account Fermi statistics and the exchange energy. The former stablilize the 
bipolaron while the latter destabilizes it. The competition of these two effects depends on 
the polaron concentration. Out treatment of the bipolaron is variational and therefore it is 
not to be excluded that some of the relative minima found here for the bipolaron energy as 
a function of variational parameters are artefacts of the method. 
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FIG. 1. The metastable bipolaron energy vs. the Coulomb coupling constant at a = 9 (in 3D) 
and a = 3.82 (in 2D). 
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FIG. 2. "Phase" diagram in the plane of the coupling constants (u 
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FIG. 3. The limiting value of the Coulomb coupling constant u vs. the polaron concentration 
for a = 7. The 3D metastable bipolaron is stabilized by the polaron environment in the region 
below the shown curve. The destabilization effect is negligible in this case. 
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FIG. 4. The limiting value of the Coulomb coupling constant u vs. the polaron surface con- 
centration for a = 3.82. The 2D metastable bipolaron is stabilized by the polaron environment in 
the region below the shown curve at a > ad ~ 43ao- The destabilization of the bipolaron happens 
for the concentrations < a < 
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